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ABSTRACT 

We present a general relativistic (GR) model of jet variability in active galactic nuclei due 
to orbiting blobs in helical motion along a funnel or cone shaped magnetic surface anchored to 
the accretion disk near the black hole. Considering a radiation pressure driven flow in the inner 
region, we find that it stabilizes the flow, yielding Lorentz factors ranging between 1.1 and 7 at 
small radii for reasonable initial conditions. Assuming these as inputs, simulated light curves 
(LCs) for the funnel model include Doppler and gravitational shifts, aberration, light bending, 
and time delay. These LCs are studied for quasi-periodic oscillations (QPOs) and the power 
spectral density (PSD) shape and yield an increased amplitude (~ 12 %); a beamed portion 
and a systematic phase shift with respect to that from a previous special relativistic model. The 
results strongly justify implementing a realistic magnetic surface geometry in Schwarzschild 
geometry to describe effects on emission from orbital features in the jet close to the horizon 
radius. A power law shaped PSD with a typical slope of —2 and QPOs with timescales in 
the range of (1.37 — 130.7) days consistent with optical variability in Blazars, emerges from 
the simulations for black hole masses M . = (0.5 — 5) x 1O 8 M 0 and initial Lorentz factors 
ljet,i = 2—10. The models presented here can be applied to explain radio, optical, and X-ray 
variability from a range of jetted sources including active galactic nuclei, X-ray binaries and 
neutron stars. 


1. Introduction 

Bulk flowing plasma close to the innermost stable circular orbit (ISCO) of a black hole consisting of 
the accretion disk and developing jet is strongly influenced by magnetic field structures, differential rotation 
of the disk, and radiative structuring. The flow velocity field and magnetic fields are dominated by the 
azimuthal and possibly the radial components. The poloidal components are non-zero owing to vertical 
structuring, random motions, radiation pressure, and dynamo action in the disk. The magnetic field strength 
generated by dynamo action in the inner disk is expected to be in equipartition with the gas pressure (e.g. 

T996] >. 


Brandenburg et al. 1995), allowing for a reasonably well developed vertically structuring in the field (Spruit 
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Variability is often observed in emission from active galactic nuclei (AGNs) with strong jet components 
such as blazars, inferred from the domination of synchrotron processes in radio to optical light curves (LCs) 
and spectra from regions at parsec scales as well as close to the base of the jet (e.g. Wagner et al. 1996; Gupta 
et al.|2012 ). Theoretical models are often applicable to emission from regions which are some distance away 


from the central region where structures are resolvable. Some of these include Doppler beaming of a stream 


or blobs of plasma accelerated to relativistic velocities along helical paths (Camenzind & Krockenbergerj 
1992 Steffen et al. [1995; Rieger 2004 1 as was recently observed in the blazar BL Lacertae (Marscher 


et al.[|2008 1, and shocks propagating along the relativistic jet (Marscher & Gea r]) 1985 1 which explains the 


variability in some blazars (e.g. |Qian et al. 1 199 4 Hagen-Thom et al.|2008| ). Quasi-periodic variability can 
be caused due to the orbital motion of the flow along helical trajectories which get beamed when their local 
angle is close to the angle to the observer’s line of sight in the inner region where the jet is just developing. 
This effect is expected to last for a few cycles. 

In a study of X-ray binaries, the inner jet appeal's to be intrinsically linked to the corona as spectral 
characteristics of its emissions are the same as those from a Comptonized corona ( Markoff et al.j2005 ). In 
another study, VLBI observations at 1.3 mm resolve structures and study the emission from the inner jet of 
M87 (Doeleman et al. |2012 1 . In the study, a direct scaling relation between the angular size and the distance 
to the object indicates that the region of emission is very compact, even inside of the ISCO (5.5 ± 0.4 
R s , where R s is the Schwarzschild radius), implying that the emitting source is on prograde orbits. These 
and similar studies indicate that there is a strong disk-jet connection, which is also likely as perturbations 
produced in the disk can be advected into the jet and amplified there via Doppler boosting (e.g. Wiita|2006]). 


Emission in AGNs such as blazars is observed in a wide variety of wavelengths ranging from radio to 
optical to Gamma-rays (e.g. Abdo et al.|2011 1 and can be inferred in many cases to arise from the jet. Intra¬ 
day variability in the optical bands is frequently observed in strongly jetted radio loud AGNs. The measured 


magnitude often changes by an order or more in less than a day (e.g. Gupta et al. 2008). Jet variability in the 
X-rays is observed on timescales of a few hours. A possible 4.6 hr QPO in a X-ray light curve of the blazar 
PKS 2155-304 and the interpretation is discussed in terms of instabilities in the disk being advected into the 
jet and intrinsic jet based processes including shocks in jets (Lachowicz et al. |2009) l. IDV timescales ranging 
between 15.7 and 46.8 ks have been detected in eight light curves from 1ES 1426+428 and PKS 2155-304 


(Gaur et al. 2010). Intra-night variability in the optical B and R bands have detected timescales ranging from 
several hours to a week in the radio loud narrow line Seyfert 1 galaxy SDSS J094857.3+002225, interpreted 
in terms of physical processes in a relativistic jet (Liu et al. ||2010 i. These and other studies indicate that 
short-term (few hours) to long-term (a few days to weeks) variability is commonly observed in the radio, 
optical, and X-ray wavelengths. A helical kinematical jet model was applied to explain the trajectories of 
blobs in the inner jet (Steffen et al7 |[ 1995 1. In the study, the model is applied to quasar 3C 345 to infer an 
inclination angle of 6.8° and typical Lorentz factors of 5.8 and 4.6 for two radio components. 


In our paper we define “blobs” to be density inhomogeneities which are acted upon by the centrifugal 
force, radiation pressure and drag, and the gravitational force of the black hole. We however do not include 
any thermodynamic evolution of the blob. The blob flows along with the bulk plasma in the region with 
a dominant magnetic field, where the plasma is constrained by flux freezing to flow along the magnetic 
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surface encompassing the jet. The variability could be attributed to orbital blobs propagating along the jet, 
the emission from which is beamed along the observer line of sight causing rapid changes in the received 
flux. Our model assumes a flow which has transited from the accretion disk onto the jet through the magnetic 
field surface anchored onto the disk at foot points close to the black hole. Thus, blobs such as orbiting spots 
could be present in the jet at various scales in the mass loaded jet. This is expected to occur regardless of 
the jet models such as the purely electro-magnetic jets or the magneto-hydrodynamic jets. The developed 
model is used to study the quasi-periodic oscillations (QPOs) phenomenon in addition to the timing study of 
simulated light curves using the Fourier periodogram and wavelet analysis. Our model can be used to study 
both the timing properties and the variable emission from orbiting blobs as well as model the jet structure 
and place constraints on the trajectories of these orbiting blobs which can be compared with observations. 


Earlier works in this direction have mainly addressed only accretion disk based variability and its timing 
properties. Models include the variability due to orbiting blobs and are applicable to optical/UV and X-ray 
wavelengths (e.g., Zhang & Bao 1991} [Mangalam & Wiita 1993). Recent models also include the effect 
of aberration based on the observer in a local static frame (Pechacek et al. 2005 >, hinting at the necessity 
for a full general relativistic treatment in the disk-jet case. A model for the QPO and its timing properties 
such as the quality factor, the break frequency inferred from the PSD shape for a general relativistic thin 
disk was presented in (Mohan & Mangalam |20l4 i. There is a necessity for a treatment of the emission 
source kinematics and light ray paths in curved space-time as the inner jet is close to the black hole; also 
general relativistic effects on this emission have not been treated in earlier works. A preliminary study was 
presented in ( Mohan|2014[ Mangalam & Mo han|2014 i. 

A schematic indicating the basic features of the wind in the context of the black hole system is shown 
in Fig. [I] The geometry consists of three zones. In the region (R s < r < Risco ) marked as Zone 1, the 
radiation pressure is dominant and drives the outflowing blobs which are sourced from a hot corona (thermal 
temperature of ~ 85 keV, e.g. Mohan & Mangalamj[ 2014j ) with zero angular momentum (as the inflow 
towards the black hole would be on plunging orbits). We consider the effects of radiation, namely pressure 
and drag in a two dimensional model in (|2j thereby extending the one-dimensional model of Abramowicz 


et al.|(1990l. In <|3j we construct a model in Schwarzschild geometry for the kinematics and emission of 


the blob for the region ( Ra < r < Rl) marked as Zone 3 and beyond (r > RjJ. The transition region 
(.Risco < r < Ra) marked as Zone 2 requires a treatment of the relativistic Grad-Shafranov to self- 
consistently describe the geometry of the field and kinematics of the flow. While we postpone this for future 
studies, in this paper we take the results from jj2]to provide the launch parameters for calculating the blob 
trajectory in a conical or funnel geometry for the magnetized jet in <|4| We assume that the blob which 
is accelerated by radiation reaches the magnetic surface where it is centrifugally driven to higher Lorentz 
factors by the co-rotation with the foot point. In the analysis in l|4] applicable to Zone 3, we calculate typical 
Lorentz factors which can be obtained both from the initial acceleration due to radiation in Zone 1 and co¬ 
rotation in Zone 2; in Zone 3 we assume that the blob has reached a final angular momentum achieved at 
the Alfven point; j x = w \D where zu,\ is the cylindrical radius of the Alfven point and Q is the spin of the 
magnetic surface. Based on the kinematics of the model, the emission is also calculated for the centrifugally 
driven flow on an assumed conical or funnel shaped geometry for the magnetic surface. While this is a 
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simplified approach which improves upon previous work, we plan to construct a fully self consistent model 
for Zone 2 in the future. We calculate the trajectory of the orbiting blob and the expected light curve from a 
special relativistic cone model as presented in Camenzind & Krockenbcrgerj ( jT992) and then from a general 
relativistic cone model and a realistic funnel model given the instantaneous position and velocity in each 
case. Results from the analysis of simulated LCs include addressing the QPO phenomenon, its evolution 
and the shape of the power spectral density for multiple emitting regions are presented in )4.4| We then 
discuss the advantages of our general relativistic funnel model and interpret the results of the simulations in 


2. Radiation driven wind in the launching region 


The effect of radiation pressure on radially outflowing relativistic particles have been studied by Abramow- 


icz et al. (19901. While Vokrouhlicky & Karas ( jl991 1 included the black hole spin on the radiation pressure 


driven radial motion of particles, the effects of thermodynamic evolution was included by Horak & Karas 


(2006). 


In addition to the effects of radiation pressure considered in (Abramowicz et al. 1990), we extend 
their one dimensional model to the general case of three dimensions (that simplifies to two for the resulting 
planar trajectories in a spherically symmetric metric) by including effects of radiation drag arising from 
the azimuthal and the 9 components of the equations of motion. The four acceleration of the blobs in 
Schwarzschild geometry is set equal to the radiation forces of pressure and drag on the blob and hence the 
Lorentz factors are derived. 


The observer geometry and the vectors is presented in Fig. 3 of (Abramowicz et aL]|1990 ). First, we 
treat the blob motion in full (r, 9, (f>) spherical coordinates. The evolution of the velocity due to the action of 
the radiation pressure and drag self consistently governs the trajectory of the blobs until it reaches the launch 
positions on the magnetic surface at larger distances. The region where this occurs is radiation dominated 
and consists of a zero angular momentum hot flow sourced from a corona. The kinematics of this region, 
Zone 1, is of interest in the current section. 

The model is cast in a Schwarzschild space-time with the line element, 


ds 2 = —(1 — 2 M/r)c 2 dt 2 + 


dr 2 


+ r 2 (d0 2 + sin 2 6d(j) 2 ), 


( 1 ) 


(1 - 2 M/r) 

where (r, 9, cf >) are the spherical polar coordinates, M = GM./c 2 is the gravitational radius with M, being 
the black hole mass. The covariant components of this diagonal metric, expressed in these coordinates are 
given by 


9a/3 = {9tt,9rr,gee,9u) = ^-(1 - 2 M/r), r 2 , r 2 sin 2 6> 


( 2 ) 
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Fig. 1.— A schematic indicating the kinematically different zones in the trajectory of a typical blob. The 
black hole with Schwarzschild radius R s = 2GM,/c 2 is at the centre. In the region ( R s < r < Risco ) 
marked Zone 1, radiation pressure is dominant and drives the outflowing blobs sourced from a zero angular 
momentum corona. In the transition region ( Risco < f < Ra) marked Zone 2, the kinematics of the flow 
are driven by radiation as well as the co-rotating magnetic field lines that are anchored to the disk. In Zone 
3, at radii Ra < r < Rl, where Ra is the Alfven radius and Ilf j is the radius of the light cylinder, the 
inertia dominates and the orbital angular momentum of the blob has reached an asymptotic value. See the 
online article for the color version. 

By writing the line element in the form 

-2 (dr 
V dt 


#\ 2 
dt J 


— r 2 sin 2 6(1 — 2 M/r) 


-i 



- r yi-2 tf/r)-(f 


(3) 


ds 2 = -c 2 dt 2 ( 1 - 2 M/r) 


1 - (1 - 2 M/r) 
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we can write the blob three-velocity components as 

1 

ft = (1-2 MM-'-* (4) 

ft = (1 - 2 

c at 
c at 

With the above representation, the proper time dr 2 = — ds 2 can be written as 

dr 2 = c 2 dt 2 7 -2 (l — 2 M/r), (5) 

where 7 = (1 — /3 2 — — /3 |) _1,/2 is the Lorentz factor. The components of the four-velocity of the emitting 

blob along the jet u a = dx a /dr = (dx a / dt) (dt / dr) are then 

u a = (-(1 - 2M/r ) _ 1 / 2 7 , (1 - 2M/r) 1 / 2 /Vy, ~Pol, —^—eP^t) > ( 6 ) 

\c r r sm 6 1 J 


up 


The covariant components of the four-velocity can be obtained by lowering the indices using the metric 

= u a g a p, 

up = (-(1 - 2 M/r) 1 / 2 7 , (1 - 2 M/r)~ 1 // 2 / 3 r 7 , r/? 07 , r sin 0(3^). (7) 


The contra-variant components of the acceleration are given by 

(jqj ^ 

+ T %,u»u v . 


The radial component is then given by 


du r du r 0 du 

a= w u+ ~ae 

du r 


u- + —u* + r^uW 


dcf) 


= ^ + —u" + —u* + T r tt ( u y + r; r ( u y + 


n m 

Using the four-velocity components from eqn. (61 and Y r tt = — (1 — 2 M/r), 
T r rr = ^(1 - 2 M/r)- 1 , V r H = -r sin 2 0(1 - 2 M/r), 


(B) 


(9) 


r 2 

a = 7 


^2 + P r (l - 2M/r)(l + PW)~^r + /5r7 2 (! - 2M/r) [Pe~^r + ( 10 ) 


A/f 

__(l_2M/r)(/3 2 + /3 2 ) + 
/3^(1 — 2M/r ) 1 / 2 


+ 


r sin# 


A,(l - 2 M/r ) 1 / 2 


3/7 


3/7 


(l + ^rV)^T + ^ Au£r + 


30 


(1 + Prl )~qq + Prl yPo~QQ + 

df3 e . 0 dy 


dfie o dPcj> 


80 


dcf) 


<90 
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The quantity mc 2 a r is the radial force where m is the mass of the orbiting blob. The 0-component of 
the acceleration is given by 


du e „ du e n du e 


a ° ~ Ik 11 " + ~m U<> + + 


dp 

- + w“ 9 + % u * + r ”*(»T + 2r -' ,r “ e ' 


<ii) 


Using the four-velocity components from eqn. (roll, T^ = — cos 0 sin 0 and T e rd = 1/r, 


a° = l- 

r 


Pr(l ~ 2M/r) 1/2 |(1 + Ph 2 )^r + Pel 2 (Pr~^r + P</>- 


dP<p\ . Pe 


dr 


Pe \ M , o2„2\ d Po , _,2 5/?r , Q dp 0 

86 


+ -\{l + Polf )-90+P<n {^-QQ + 


. r , n , a2 2^Pe , a 2 f a d Pr , 0 dP* 

+ V^~e\ {l + M) ~dp +M Y r ~dp + ^~dp 


cot 0 


H- 

r 


( 12 ) 


The quantity mc 2 ra e is the force in the 0 direction. The azimuthal component is given by 


cr = 


dvP 

dr 

duP 


u r + 


dvP a du$ J. ^ 

W u+ w u+r > u 


lk u + 1 w u ‘ + + 2 r 7" r “* + 


(13) 


Using the four-velocity components from eqn. ( 61 , T'^7 = 1/r and T^V = cot 0, 


00 


or = 


7 


r sin 0 


Ad - 2 M/r)^ {d + /3j7 2 )^ + A ~ + 7 2 (ft-f 2 + A^)) } 


+y {(! + Ad// + 2 A A7 2 + a//) J 


+ 


r sin 0 


■2„,2\ d P<t> , a „,2 (o d Pr , 0 9/3 e 


(1 + } ^ + /V7 ( 


(14) 


The quantity mc 2 raP represents the force in azimuthal direction. Radiation force drives the accelerated 
outward motion of the blob while the initial azimuthal and angular motion however is continuously retarded 
by the effect of radiation drag of an assumed spherically symmetric radiation field. 


Following Abramowicz et al. (1990), the contra-variant components of the radiation energy flux are 
given by, 

F a = h*T^u u , 


( 15 ) 



























where /i“ = 5" + u n u /l is a projection tensor and T ,l,J is the radiation field energy density. The tetrad 
components of the energy density tensor can be evaluated using 


T (0(fc) = / /( r ) n «n (fe) dU 


(16) 


where I (r) is the radiation field intensity, are the unit vectors describing the photon trajectory and dil 


is an element of the solid angle subtended on the sky of the local observer (see Fig. 3 of Abramowicz et al. 

(17) 


1990). The energy density can thus be evaluated from this using 

T tt = _ 2 M/r)- 1 = 27r/(r)(l - 2M/r) _1 (l - cost/) 

T tr = T rt = T (t)(r) = T (r)(t) = 7r /( r ) sin 2 ^ 


T rr = rWW( 1 - 2 M/r) = 2 ^l(l - 2M/r)(l - cos 3 rj) 
eo_TVW) _irl(r),_ 3 


rj-iUt7 _ 


rjHjxj) _ 


3 r 2 


(cos 3 r/ — 3 cos rj + 2) 




7r/(r) 


- (cos 3 77 — 3 cos r/ + 2), 


r 2 sin 2 0 3r 2 sin 2 0 

where 7? = 6M is the radius of the emitting surface, rj is the viewing angle from the zenith for the stationary 
observer and is related to 5, the zenithal angle at the point of emission, R by 

M (1 - 2M/r) 1 / 2 


snip = 


where sin 5 = 1 for R/M >3/2 and sin 6 = 


r (1 — 2M/7?) 1 / 2 
3\/3 M 


sin 5, 


(18) 


1990). 


2 R 


(1 — 2 M/R) 1 / 2 for R/M < 3/2 (Abramowicz et al. 


If a is the cross section of the orbiting blob over which the radiation force acts, the quantity aF a 
represents the radiation force. The radial component of the energy flux is given by 

F r = [(1 + u r u r ){T rt u t + T rr u r ) + u r u t (T tt u t + T tr u r ) + u r u e (T ee u e ) + u r u^T^ U( p)\. (19) 

Using the components of u a from eqn. ([6]), u a from eqn. ([ 7 ]) and T' w from eqn. ( [ 17 ] ), the radiation 
force in the radial direction is given by 


aF r = 7r/(r)cJ7 3 (l — 2 M/r) 1//2 


1 + ^2 ^ Sin2 7 + \PrO- - COS 3 r/) 


( 20 ) 


-/3 r ((3r sin 2 p - 2(1 - cost?)) + ^{0e + /3|)(cos 3 rj - 3cosp + 2) 


The 0-component of the energy flux is given by 

F 9 = [(1 + u 9 uq ) T 99 uq + u e ut{T tt u t + T tr u r ) + u e u r (T rt u t + T rr u r ) + 


( 21 ) 
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Using the components of u a from eqn. (|6j), u a from eqn. ([7]) and the components of T ,w from eqn. 
([17]), the 0 -component of the radiation force is given by 


aF e = 7r7(r)cr 


1 + /3gT 2 

37 2 


(cos 3 T] — 3 cos ?7 + 2) 


( 22 ) 


+2(1 — cos 77 ) — -(3 2 ( 1 — cos 3 ?/) + /3|(cos 3 ?/ — 3 cos 77 + 2) 


The azimuthal component of the energy flux is given by 

= [(1 + + u*u t (T tt u t + T tr u r ) + u^u r (T rt u t + T rr u r ) + u*u e T 6e u e \. (23) 


Using the components of u a from eqn. (|6j), u a from eqn. ([7]) and the components of T ,IV from eqn. 
© the azimuthal component of the radiation force is given by 


aF * - 7 r/ ( r ) g 7 3 / 

r sin 0 


+ /%7 2 

37 s 


(cos 3 ?7 — 3 cos 77 + 2 ) 


+ 2(1 — cos?]) — ~/3r( 1 — cos 3 r/) + /3|(cos 3 ?/ — 3cos rj + 2 ) 
o 


(24) 


2.1. Dynamics based on radiation pressure and drag 


The Eddington parameter is the ratio of the disk luminosity L to the Eddington luminosity LEdd and 
is given by, 

L CTTl(R)aR 2 (l - 2 M/R) 1 / 2 


T f = 


7>Edd 


GM,mc 


(25) 


where +Edd is corrected for the general relativistic redshift factor (1 — 2 M/R) 1 / 2 and a is the cross 
section of the orbiting blob over which the radiation force acts and 


7(0 

HR) 


1 - 2 M/R\ 2 
1 - 2 M/r ) 


(26) 


is obtained from the energy conservation along the null trajectory, in the stationary frame (Abramowicz et al. 


1990). 


The first equation of motion is from the radial components of the acceleration and the force. The radial 
component of the force imparts the radial acceleration onto the orbiting blob, i.e. aF r = mc 2 a r . Using 
eqns. © and ( [20] ) and the representation x = r/M and X = R/M, we can eliminate the dependence of 
the equation on M. Then, expressing the equation in terms of the Eddington parameter T f using eqns. 


and ( 261 , 
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5/3, 


M2 + A-(! - 2 /^)( 1 + Prl )-Q^ + Prl (1 - 2 /®) \ Po-fa + 

dpr 

d6 


d/3 e d/3<f> 


dx 


- 1(1 - 2/x)(ft + ft) + w y \ (1 + n W) -XL + A7 2 


A)(l - 2/x) 1 / 2 




+ 


Mi 2/x) 1/2 ; (i + ftV) ^. + ^ 2 


x sin (7 

r e7 /i - 2 /x\ 3/2 
^ 2 " ^ 1-2/sy 


dPo 5/fy\ 


1 yy 7 (2 sin 2 J, + lft( 1 - cos 3 17) 


/ «2 + fl2 > 

— / 3 r ( f 3 r sin 2 77 — 2(1 — cos 7 ) H— 6 -— —-(cos 3 ?/ — 3 cos 77 + 2 ) 


(27) 


d/ 3 g_ 

’ 50 


’ 50 / 


The second equation of motion is from the 0-components of the acceleration and the force. The 0- 
component of the force imparts a drag on the acceleration of the orbiting blob and aF e = mc 2 a e . Using 
7r = (l - ft 2 ) 1 / 2 , eqns. ([12]) and \2l\ and the representation x = r/M and X = R/M and expressing the 


equation in terms of the Eddington parameter T e using eqns. ( [23] ) and 

MI- 2/x) 1/2 {(1 + AV)§| + At 2 (F// + A(2’) + * 


+ ^Ui + 3h 2 )^+M 2 U^/ + , 


X 


dp# 

’He 


+ 


,2 2^Pe , 0 , 2(0 df3r 


xsind | (1 + ^V)^ + M 2 (&^+, 


dP<j, 

’ 5y> 


cot 0 


r e7 A> 


1 - 2/X 


X 2 (l - 2/Xy/ 2 \l-2/x 

|ft(l-COS 3 ft + l(J '002 


2(1 - cos?/) - 7 / 3 2 (l - cos 6 rf) + 7: ( + 2 / 3 ^ ) (cos 3 rj - 3 cos?/ + 2) 


(28) 


The third equation of motion is from the azimuthal components of the acceleration and the force. The 
azimuthal component of the force imparts a drag on the azimuthal acceleration of the orbiting blob and 


aF° = rnc 2 (F. Using eqns. (14 1 and (24) and the representation x = r/M and X = R/M and expressing 

























-li¬ 


the equation in terms of the Eddington parameter T e using eqns. (25) and 


/5 r .(l-2/x) 1/2 |(l + /3|7 2 )^+/3 0 ^-cot0^+ 7 2 ^/3 r .^ + /3^ ) ) } (29) 

+ §{ (1 + Pll 2 )- + W* cot 9 + /^ 7 2 + Po^f) } 




+ 


. o2 2\^B(j> , a 2/^o &Br n dfie\ 


x sin 0 

_ T e7 ^ / 1 - 2 /X V 
X 2 (l - 2 /X ) 1 / 2 V 1 -2/x ) 

2(1 - cos 7 ) - |/3 2 (1 - cos 3 7 ) + \{\+ 2/3g ) (cos 3 7 - 3 cos 77 + 2 ) 
3 3 \7r 


2.1.1. Purely radial motion 


In the limit B<f> = 0 and Be = 0 , i.e. purely radial outflowing motion, only the first equation of motion 
is relevant. The equation reduces to 


(1 - 2/xHp r d A + 4 = 4 

ax x z x z 


1 - 2/X \ 3/2 

1 - 2/x J 


— sin 2 7(1 + /3 2 ) + -fl r ( 4 — cos 3 7 — 3 cos 7 ) 

O 


, (30) 


which is the same as that derived in (Abramowicz et al. 1990). Solving for Br and hence y r = (1 — 6/) */ 2 
as a function of the initial launch velocity, Br,i and x indicates a constant saturation value of 2 — 7 for 
T e = 0.3, PrA = 0.86 — 0.99 at distances r beyond ~ a few tens of M from the black hole and is plotted in 
Fig. [2] In the simulations, changes to T e did not change these results. 

For the range of Br,i = 0.86 — 0.99 and x = 6 — 20 , the 7 r contours saturate at constant values in the 
range 2 — 7 for large x, consistent with the simulations of Abramowicz et al. (1990). 


2.1.2. Poloidal motion 


A poloidal outflow in the radiation region (B<f> = 0 , do B 0 , Br B 0 ) is trajectorially equivalent to the 
case of motion on the equatorial plane (6 = ir/2,Be = 0 , Bcj> B 0 , B r B 0 ) due to spherical symmetry 
of the Schwarzschild geometry. We appeal to the zero angular momentum nature of the outflow source in 
the corona region (see Fig. [TJ to motivate this poloidal flow. Thus, we are left with two coupled partial 


differential equations to solve, modified versions of eqns. (27 1 and ( 281 with (B<j> = 9- Bo B 0 , Br B 0)- The 
bulk Forentz factor 7 = (1 — Br ~ Be ) 1//2 in this case and j r = (1 — /3 2 ) -1 / 2 ; 70 = (1 — ^) -1 / 2 . The 
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Fig. 2.— Contours of 7 r as a function of the initial launch velocity P rA and x. The x r value tapers off to 


constant values at large x for = 0.86 — 0.99 consistent with the simulations performed in Abramowicz 


et al. (1990). The simulated x r are in the range 2 — 7. 


resulting equations are, 


1 


— + Prl (1 - 2 /x) 


x \7e 90 


1 9 fir , a a 9 ' 6 0 

— 2 7 )- ^ PrPe-X— 

7 g ox ox 

9 Pe 

' do 


-^(l-2/X 


(31) 


r e7 (i - 2/x \ 3/2 r i 


FU-! /xj 

R 2 

— /3 r ( /3 r sin 2 ?7 — 2(1 — cos 77 ) + (cos 3 77 — 3 cos 7 + 2 ) 

O 


for the radial motion and the 0 -equation is, 

■' 2 aa - 2 d/3 r /3 e \ Pel 2 ( 1 9 Pe 


PrO- - 2 /x ) 1/2 ^ + PePrl 2 


dx x / x Vt 2 dO 




TelPe 


X 2 (l — 2/X) 1 / 2 \ 1 - 2 /x J 


1-2/lV 


(32) 


2 1 

2(1 — cosr/) - p 2 (l — cos 3 7 ) H-^(cos 3 7 — 3 cos 7 + 2 ) 

3 37 / 


In the simulations carried out, we set F f = 0.3. The initial value p rA was varied between 0.01 - 0.99, 
the range being chosen based on the P r range simulated in (Mohan & Mangalam j2014 ) for a relativistic thin 
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disk ( Shakura & Sunyaev|1973 ; Novikov & Thome|1973 ) in the context of the quality factor Q observable 
in the inner region of the accretion disk very close to the innermost stable circular orbit. There, we obtained 
the range 1.36 x 10 5 — 0.99 for a range of the disk viscosity parameter a = 0.01 — 0.4 for r > 6 M — 20 M. 
The range of initial velocity /3g ti is chosen based on the argument that if the motion is along Keplerian 

orbits, Po = ( 1 — 2 M /r )^ 1 ' 2 - — = (1 — 2 M/r)~ 1 ' 2 -£l(r). As Q(r) = c(Af/r 3 ) 1 / 2 for Keplerian orbits, 
fjg = (1 — 2M/r) _1 / 2 (Af/r) 1 / 2 . We thus evaluate Pg i to be in the range 0.25 — 0.61 for r = 6 M — 20 M. 


The initial values p r ^ and Bq a were taken in the range 0.86 — 0.99 (similar to the case of purely radial 
flow) and 0.25 — 0.61 respectively. The initial values of the variables x and 6 were taken in the range 6 — 20 
and 0.01 — 7t/2 respectively. The contours of fig and 7 = (1 — /3 2 ) -1 / 2 where P = (/3 2 + Pg) 1 ^ 2 as a 
function of x and 6 are plotted in Fig. [3] Simulated final values of Bo are in the range 0.18 — 0.81 and 7 are 
in the range 1.1 — 26.3, higher than in the purely radial case though, 7 saturates to low values for smaller 
initial values of x. It is inferred from the simulations that Pg and 7 decrease rapidly; this decrease is large 
for smaller initial values of x, while their decrease is more gradual at larger initial values of x. These effects 
can be attributed to the drag force due to the azimuthal component of the radiation pressure which plays an 
important role in stabilizing the saturation values of 7 at an early stage even if the source of the outflow had 
an initial angular momentum. The simulations thus indicate that the flow becomes rapidly radial. 


We plan to present details of classification of flows and hence, a more comprehensive exploration of the 
parameter space and the resulting trajectories in a paper in preparation. Here, we have derived the typical 
values of the poloidal P p and 77 that are the inputs to typical initial values for Zone 3 ( Ii/\ < r < R[J and 
beyond (r > Rl) where the angular momentum is set by its value at the Alfven point. 


3. Trajectory and emission geometry in and beyond Zone 3 

From the metric in eqn. as there is no explicit dependence in the equations of motion on the cf) 
and t coordinates, there are two Killing vectors associated with this geometry given by £ = (1,0,0,0) and 
77 = (0, 0, 0,1). If u is the four-velocity of the emitting blob, the Killing vectors can be used to evaluate the 
constants of motion for the Schwarzschild metric including the the total energy 

u t 

e = - < " 1= (1-2M/R)' <M) 

and the conserved angular momentum 

j = —rj ■ u = VLR 2 sin 2 9, (34) 


where u l is the time component of the four velocity u and fl is the angular velocity of Keplerian orbits 
in Schwarzschild geometry. The condition u • u = — 1 in case of the Schwarzschild metric gives 


U* = 7 jet = 1 - 2 M/R - 


R 2 /c 2 


-1/2 


- R 2 6 2 /c 2 - R 2 sin 2 0 (j) 2 /c 


12 


1 - 2 M/R 


(35) 
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Pe (x,0) 



x = r/M 


Y(x.0> 



x = r/M 


Fig. 3. — Contour plots of pg(x,9) and 7 (x,9) = (1 — /3 2 ) -1 / 2 where /3 = (/3 2 + Pg) 1 ^ 2 for P r ,i = 
0.86 — 0.99 and = 0.25 — 0.61. Left plot: simulated /3 g are in the range 0.18 — 0.81. The decrease in 
Pe is large for smaller initial values of x. Right plot: the decrease in 7 also follows the same expected trend 
to that of fig. The simulated 7 are in the range 1.1 — 26.3 higher than in the purely radial case though with 
saturation of 7 to low values occurring even at smaller x. These effects can be attributed to the drag force 
acting to rapidly cause the loss of any initial angular momentum as well as stabilize the flow at small x. 


This can be used to write the general expression for the velocity Pj et of the emitting spot as 

R 2 /c 2 


Pjet — 


— R 2 9 2 /c 2 — R 2 sin 2 9 <p 2 /c 


2 /„ 2 


1 - 2 M/R 
= (l — 2M/R — 1/ "fjet ) 1/2 


(36) 


The various emission and direction vectors include the radial vector n, the direction vector of the 
initially emitted light ray ko and the final direction vector pointing along the observer line of sight k. A 
general path showing the source motion along a helical trajectory and the emission geometry along with the 
above vectors is presented in a schematic in Fig. [4] 

The vector ko is given in terms of k and n as 


, sin a, sin(-ib — a) 

k 0 = -—rk H -;—-—-n; 

sin y) sin yi 


(37) 


this identity can be verified applying dot and cross products of k and n and using their geometry given in 

Fig- 0 
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Fig. 4.— General helical flow geometry showing the path of the emitting source on the jet (blue) launched 
from the cylindrical radius wq emitting a light ray which is subjected to the light bending effect (red). If 
the time of emission is t em , the time at which the signal is observed t 0 b s = {t em ~ tadv + At lb) (1 + z) 
where t a d v is the time taken for the source to traverse the distance R cos t/j and At lb is the correction due 
to light bending effect in curved space-time. The blob trajectory is shown in blue, emission vectors in black 
and geometrical quantities including distances in grey. See the online article for the color version. 
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In the Cartesian ( x , y, z) coordinate system centred on the black hole, the components of the direction 
vector k are given by 


k = (sini, 0, cos i) (38) 

and the components of the emission vector n are given by 

n = (sin 9 cos </>, sin 9 sin <p, cos 9 ). (39) 

The initial emission angle a can be written in terms of the initial direction vector ko and the emission 
vector n as cos a = ko • n. The final angle of the emitted ray ib can be written in terms of the final direction 
vector k and the emission vector n as cos ip = k • n. This can be expanded as 

cos ip = cos i cos 9 + sin i sin 9 cos cp. (40) 


Assume that an emitter located at a radial distance R emits a light ray at an angle a with respect to 
the radial vector. The light ray undergoes bending due to space-time curvature to emerge at a final angle 
-ip with respect to the radial vector and the relationship between ip and a can be determined by solving and 
combining 

ip = —r(l/6 2 — l/r 2 (l — 2M/r))^ 1 / 2 (41) 

Jr 

and 

sina = 4(1-2M/ii) 1 / 2 , (42) 

IX 

where b is the impact parameter shown in Fig. [4j The derivation of the above equation is presented in 
Appendix [A] This propagation effect can be approximated by the expression 

1 — cosck = (1 — cosy))(l — 2 M/R). (43) 


([Beloborodov|2002) 


Using the above approximation, we find 

= (1 - 2M/R) 1 / 2 (1 - 2M/R+{4M/R)(l + cosip)- 1 ) 1/2 . (44) 

sm^i 


The travel time will be different for photons emitted from different radial locations R. The difference 
between the travel time for a ray emitted from a position R and that from a position at the centre of the 
coordinate system can be obtained from the integral 


At 


f 00 dr 
r 1 - 2 M/r 



r (l - 2 M/r) 



( 45 ) 


The time delay due to light bending can be approximated (Poutanen & Beloborodov|20061 as 
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A^lb = 



(46) 


where y = (1 — cosip) and 6 = 2 M/R. The time delay is defined with respect to straight paths from the 
point P to the observer as shown in Fig. [4j 


The time of advance due to orbital motion of the emitting source, t a d v can be expressed in terms of 
geometric factors from Fig. [4] as 


tadv = (1 - y)R/c- (no • \)wq/c, (47) 

where the constant shift (no • i) accounts for the initial offset in R in the projected i direction towards the 
observer. If the time of emission in the source frame is t em , the observed period of the signal t Q b s is reduced 
due to the effects of the component of the source moving towards the observer line of sight and due to the 
time delay induced by light bending and the disk inclination and is given by 

tobs = [tem — tadv + A^bb] (1 + z) — [tem + (y — l)R/c + (no ■ ijroo/ c + A£lb] (1 + z) (48) 
= [t em ~ {R/c) cosip + (n 0 • i)o7 0 /c + At LB \ (1 + z). 


This is used to map the time of emission to the time of observation by accounting for light bending where 
At lb is the time delay from eqn. (461 due to light bending and where z is the cosmological redshift of the 
emitting source. The effect of the time delay due to the light bending and disk inclination is to change the 
phase of reception of the observed signal. The phase change is expected to systematically increase as the 
emitting source moves along the expanding jet towards the observer. Also, we set the constant offset so that 
t 0 bs(tem = 0) = 0 for no light bending. 


4. Flow trajectories in and beyond Zone 3 and resulting light curves 


We now discuss the construction of the light curve for the variability model. The motion of the flux 
frozen blob is along a magnetic surface with foot points anchored on the accretion disk. A basic schematic 
of the blob motion along the magnetic surface is presented in Fig. [5] 

This surface co-rotates along with the disk at an angular frequency given by 


n cM 1 ! 2 

F = r4 2 + aMW' 


(49) 


where a is the black hole spin parameter and wf is the cylindrical radius to the foot point from the coordinate 
system centred on the black hole (z = 0). The radius of the light cylinder ujf is the position on a cylindrical 
surface along which the plasma moves at the speed of light and calculated as 

3/2 

C ZU rp 

WL= n F = uvi + a 


(50) 
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Fig. 5.— Helical trajectory of an emitting blob in Schwarzschild geometry, constrained along rotating 
magnetic field lines with footpoints on a Keplerian disk (at cylindrical radius w 0 ). In the cone model (left 
plot), the half opening angle of the jet is 9q. In the funnel model (right plot), the flow is asymptotically 
bound by a cylinder of radius wj at large z. The blob trajectory is shown in blue, emission vector in red, 
geometrical quantities including distances in grey and the jet magnetic surface shape in green. See the online 
article for the color version. 
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The cylindrical radius from which the blob is launched is given by vjq = /vjl where / = 0.1 - 10. The 
specific angular momentum of the blob, j is conserved along the magnetic surface. This occurs beyond the 
Alfven point with a cylindrical radius voa = where x\ < 1. The conserved angular momentum is 

then given by 

3 = Joe = x 2 A zujn F . (51) 


The Hr in eqn. (49) is the angular frequency of orbits in the Ken - metric and we make use of a / 0 only 
for comparison with the special relativistic model presented in Camenzind & Krockenberger ( j!992) l. For the 
general relativistic models, we make use of the Schwarzschild geometry (a = 0) to maintain consistency 
with the treatment of the Doppler boost factor, aberration and light bending. 

The instantaneous position of the emitting spot is given by 


x s = ( x,y,z ) = cos w{t) sin (/>(£), z(t)), (52) 

where zu is the cylindrical distance to the emitting spot from z(t). In our variability model, a prescription 
for w = w{z) is given by the geometry of the magnetic surface. The velocity components (x s = dx s / dt) 
are then 

x s = ( oc , y 1 z) = (vo cos 4>(t) — , t±7 sin + w(f>(t),z). (53) 

The angle between the velocity vector of the spot x s and the initial direction vector ko is given by 

x s • ko 1 / sin a . sin(i/i — 

cos £ = .. . = tt-t -— ~x s ■ k 4- ; —- 

|x s | |x s | \siny smip 



Using the initial direction vector ko from eqn. ([37]), the final direction vector k from eqn. (38) and the 


radial vector n from eqn. (39), we can write the most general form of cos f; as 


cos£ = — 


1 ( sin a . 


x s | \sin^ 


x s • k + 


sin('i/) — a) . 


sin V’ 


x s • n 


(55) 


1 J sin ol siYiiip — ct) 

= ——r < - [w cos 4> sin i — wH sin <f> sin i + z cos i\ H- - —-i[tan 6 q sin 6 + cos 6] 

sin ip sin ip 


|x s 

1 

|x s 


sin ex. sin( , 0 — 

[tan 6q cos 4> sin i + cos i\ H- : —-—- [tan 6q sin 6 + cos 9] 


sin if) 


sin if; 


^ sin a . , . . 

—mi l -- sm (b sm * 


sin ijj 


The Doppler factor D evaluated in an instantaneous stationary frame at the source using the emission 
vectors given in eqns. (37), (38) and (39) is given by 


D = 


1 

7jet(l - Pjet COs£)‘ 


(56) 


In the observed frame, the above expression must be modified to account for the gravitational redshift. 
Thus, the effective redshift factor g which is the ratio of observed to emitted energy of the emitted ray must 
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include the above Doppler factor as well as the gravitational redshift effect and is given by 
g = = (1 - 2 M/R)'*D = {l 7, 2M J R)l/ \ v 

-^emitted 'T jet (1 Pjet COS £) 


(57) 


An alternate derivation of the above equation is presented in Appendix [A] If 7 j e t,i is the bulk Lorentz 
factor of the flow at the initial launch radius vjq, the energy conservation condition along the trajectory can 
be written as 


ut _ 7jet _ _ ljet,i _ 

(1 - 2 M/R) (1 - 2 M/R) " (1 - 2 M/w 0 ) ' jetJ 


(58) 


This can be used to express 7 j et in terms of the initial parameters wq and Jj e t,i- From the simulations 
invoking radiation pressure and drag in ( 2.1 we motivate stable values of 7 j e t,i ~ 2 — 6 at larger x such 
as the initial launching point (foot points on the disk) which we have taken to be between ( 0.1 — 10 ) light 
cylinder radii. The light curve is given by the spectral flux density observed F v (t). If F' v is the spectral flux 
density in the co-moving frame, these arc related by the expression 


F u (t) = g x (t)F u (t), 


(59) 


where A = 3 + T for a resolved blob of plasma and A = 2 + T for a continuous flow; T is the spectral 
index which is the slope in the relation F v oc u T between the spectral flux F v and the emission frequency 
in the observer frame. Either case of A is possible depending on the particular application to observations, 
but for the purposes of illustration, we take A = 2 + T in order to compare our results with Camenzind & 
Krockenberger j ( | 1992j ). 


As t is the time in the source frame, we obtain the time in the observer frame using the transformation in 
eqn. (48 1 with the time delay AtLB to obtain F u = F u (t 0 b s ) in the observer frame. The beaming effect will 
be prominent in the light curve when the angle £ between the initial emission vector and the velocity vector 
of the source is close to the disk inclination angle i between the normal drawn in the coordinate system 
centred on the black hole and the observer direction. The quasi-periodic behaviour from orbital signatures 
when beaming occurs is expected to be distinctly visible in observed light curve and could last for only a 
few cycles. 


4.1. Special relativistic cone trajectory 


We construct a special relativistic jet model based on that presented in Camenzind & Krockenberger 


(1992). The model consists of a relativistic blob in a cone geometry as shown in Fig. [5] The kinematical 
prescription includes expressions for the cylindrical radius w and the associated velocity vj in terms of z 
and i. In the cone model with a jet half-opening angle Oq, 


-CD = wq + z tan 6q. 


( 60 ) 
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Using the above equation, the velocity associated with zu(t) is given by 

zu = z tan . 


(61) 


The conservation of angular momentum gives an azimuthal velocity <f> = Q(t) = j 00 /zu 2 and the 
phase is given by (pit) = f* (l(t)dt. The condition u • u = —1 in case of the Minkowskian metric, an 
approximation used in Camenzind & Krockenberger (1992 1 then gives 


l + ui 


^ et i-w 2 n?/c 2 


l+ui 


1 - 


7 2 

J oo 

zu 2 c 2 


(62) 


where u p = u t (+ 2 + i 2 ) 1 + is the poloidal velocity. In the region where collimation of the jet occurs 


(zu » wl ), the bulk Lorentz factor 7 j et becomes a constant due to energy conservation as eqn. (33 1 can 
be written for large R as 7 j et = e. At large R, zu is also large. Hence, u p tends to a constant value and 
like (Camenzind & Krockenberger| |T992] ), we make this approximation. For a jet half-opening angle Oo, as 

i = cu p cos Oo /7 jet., we obtain 


Ur. 


Z = 


1 + u 2 


7 2 
Joe 


cos Oo c 2 - ^I 


VO* 


1/2 


(63) 


The cylindrical distance vo = zo(t) in this particular case due to the simplification of u p being a 
constant. From eqn. w = z tan Oq . Using z from eqn. (|63j) in this, 


CUr 


VO = 


1 + u 2 


: Sill ( 


,-2 X VS 
^2 _ J oq 

W 2 


(64) 


The above expression can be written in terms of w and zo as 

CUp sill Oq 


'UJ 'UJ 




1 + u 2 


(65) 


On integrating the above equation, we obtain 


zo(t) = ^ | 1 + 



2 \ 1/2 


sin Oot 


( 66 ) 


The magnitude of the source velocity is given by 

|x s | = (x s • x s ) 1/2 = (zu 2 + zo 2 Vl 2 + i 2 ) 1/2 = (z 2 sec 2 0 o + j^o/ru 2 ) 1/2 (67) 
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In the absence of the light bending effect, tp = a in which case the direction vectors ko = k. The angle 
between the velocity vector x s and the initial direction vector ko, cos£ from eqn. ( [55] ) then reduces to 

cos£ = ——r(i(tan 0 o cos^sini + cosi) — zuDsin^sinz). ( 68 ) 

Xs | 


The velocity z can be obtained from eqn. (631 and w from eqn. ( 661 . The effective redshift factor g(t) is 
then given by 


1 


9 = 


7jet(l - P jet COS £)' 


(69) 


We can use cos £ from eqn. ( 681 , 7 j et from eqn. (62), 3 ]e t = y 1 — 1 /yjy, and w from eqn. (661 in the 
above equation to obtain g = g(t ). The light curve is given by F(t) = g 2 +r (t). As t is the time in the 
source frame, we obtain the time in the observer frame using the transformation in eqn. (48) without the 
time delay At lb i n this case to obtain F = F(t 0 b s ) in the observer frame. 


To illustrate the validity of the model developed in (j3] we simulate a light curve for a general BL Lac 
object using the expressions presented in this section. For this, we consider the same parameters which were 
considered in Camenzind & Krockenberger (1992), i.e. u p = 3, * = 0 05°, M. = 5 x 1O 7 M 0 , a = 0.8 and 
wo = 10wl- The parameter u p = 3 is chosen in order that 7 j e t,i at large w ~ 3.2 which is within the range 
predicted by the effects of radiation pressure and drag in § 2 . 1 |and is in the same range as that considered in 


Camenzind & Krockenberger (1992 1 . The resulting light curve is plotted in Fig. [ 6 ] which is similar to Fig 3 


of Camenzind & Krockenberger (jl992) with phase shift. 


4.2. Fully relativistic cone model 


We construct a fully relativistic jet model in Schwarzschild geometry which consists of a relativistic 
blob in a cone geometry as shown in Fig. [5] Here we consider general relativistic effects which include time 
delay due to light bending, the Doppler and gravitational redshift and the treatment of the poloidal velocity 
as a general function of the geometrical and kinematical parameters as opposed to the constant assumed in 
the previous case. The kinematical prescription includes expressions for the cylindrical radius w and the 
associated velocity w in terms of z and z. 


The cylindrical distance w{z) is given by eqn. (60) and the velocity associated with w(z) is given by 
eqn. ( [ 6 l] ). The radial distance R(z) is given by 

R(z) = ( w 2 {z ) + z 2 ) 1 / 2 = (wq + z 2 sec 2 9q + 2 zi 7 o 2 tan#o ) 1//2 (70) 


and the velocity associated with R(z) is given by 


1 z 

R = — [ww + zz) = — (wq tan 0q + z sec 2 60 ) ■ 
R R 


(71) 


The polar angle 0(z) is given by 


sin0(^) = w{z) / R(z) 


(72) 
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Observer time t a b s (days) 

Fig. 6 .— Simulated light curve for u p = 3, xa = 0.9, i = 0.05°, M. = 5 x lO 7 Af 0 , a = 0.8, roo = 10 vol 
and T = 1. The light curve (red curve) shows the quasi-periodic oscillation expected from a general BL Lac 
object when the special relativistic formalism is used in the simulation. See the online article for the color 
version. 


and the velocity associated with 9{z) is given by 


z\ R I R 


(73) 


The conservation of angular momentum gives an azimuthal velocity f = Q(z) = joe / zu‘ 2 . Using the 


condition u • u = —1, from eqn. (351 


U = q /jet = ( (1 - 2 M/R) - r2 ^ 2 ( TUg + 


z 2 f 2 (roo tan^sec 2 9q) 2 \ J, 


(1 - 2M/R) 


i 2 

too 


VJ 2 C 2 


- 1/2 


(74) 


From the energy conservation condition along the trajectory, 


u 






(1-2 M/R) (l-2M/ro 0 )’ 


(75) 


where z jet , is the initial Lorentz factor. Using 7 j et * from eqn. (62) for the case of w = zuq at z = 0 and a 

f \V2' 


constant u p , we obtain 7 j et) i = 


1 + U r, 


1 ~ 


7 2 
Joo 

2„2 


and 


UJ q C' 


t (1 - 2M/R) 

u ~ ljet ~" /jet ’\l-2M/w 0 y 


( 76 ) 
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Using the above expression for if in eqn. (74 1 , we obtain an equation for i given by 


z = cr(( 1 - 2 M/R) ~Sh~ W2 (*4 + + 2sec2<, »> 2 ' " I/2 


tzr2 c 2 (1 - 2M/R) 2 r /j et , l 


(1 - 2 M/R) 


(77) 


In the above equation, w = zu(z) from eqn. (60 1 and R = R(z) from eqn. (70). The above expression 
can be integrated to obtain z = z(t) as 


/ 


r* dz 

lo 


2 (tuo tan 6 q + z sec 2 do) 2 
W ° + (1 - 2 M/R) 


\ 


1/2 


(1 - 2 M/R) - 


jL (1 — 2M/zuo) 2 1 


/ 

Jo 


= c dtt. 


(78) 


TO 2 C 2 (1 _ 2M/R )2 ^ . / 


Once we obtain z = z(t), this can be used to evaluate i = i(z(t)), tn = w(z(t)), w = za(z(t)). The 
magnitude of the source velocity is given by 


|x s | = (x s • x s ) 1//2 = (do 2 + zu 2 Ci 2 + i 2 ) 1//2 = (z 2 sec 2 0 O + j^/ro 2 ) 1 / 2 . 


(79) 


The velocity z = z(z(t)) and all other terms in the above expression can be cast in terms of z(t). In the 
current model, we consider the light bending effect in which case the direction vector ko is obtained from 
eqn. © and the angle between the velocity vector x s and the initial direction vector ko given by cos £ is 


given by the general expression in eqn. (55 1 . 


The effective redshift factor g(t) is then given by the general expression in eqn. (57). We can use cos £ 


from eqn. (55 l, z je t from eqn. (74), j3j et from eqn. (361 in the equation for g to obtain g = g(t). The light 
curve is given by F(t 0 b s ) = g 2+r (t 0 b s ), expressing t in terms of t 0 & s . The importance of making use of a 
fully relativistic model when compared to the earlier special relativistic model is presented in terms of the 
difference between the t 0 b s in both cases. It is seen that for the fully relativistic model, t n i )fi tends to increase 
systematically with t em . This difference in timescale is of the order of ~ 5 days by the end of the simulation. 
This difference is plotted in Fig. [7] which justifies the use of the fully relativistic model. 

To illustrate the model developed in jj3] and the comparison with the special relativistic cone case in 
§4. 1 [ we simulate a light curve for a general BL Lac object using the expressions presented in this section. 
Using again u v = 3 as an initial condition, with i = 0.05°, M, = 5 x 10 7 Mq, a = 0.8 and wq = 10 wl, 
the resulting light curve is shown in Fig. [8] The parameter u p = 3 is chosen in order that jj e t,i at large w ~ 


3.2 which is within the range predicted by the effects of radiation pressure and drag in (2.1 and is in the 


same range as that considered in Camenzind & Krockenberger (1992). This is a comparison between the 


light curves in the fully relativistic cone case and the special relativistic case which indicates an increased 
amplitude in the former; there is also a systematically increasing phase lag compared to the latter case as 
expected due to the effect of time delay caused by light bending. The ratio of the amplitudes gives a 12 % 
increase in the maximum amplitude due to the GR boost factor in g. The beaming effect is observed for the 
last two cycles for the general relativistic cone case where there is an increased amplitude which then settles 
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Emission time t em (days) 


Fig. 7.— Curves showing t Q b s versus t em for the fully relativistic cone model (blue) and the special rela¬ 
tivistic cone model (red, dashed). The t 0 b s in the former increases systematically and reaches a difference of 
~ 5 days by the end of the simulation with respect to the latter. See the online article for the color version. 


down. The increased amplitude ratio and the phase lag change can be attributed to the integrals over time 
of intensity and phase amplifying the small but early differences in phase. Even a small phase difference 
between the general relativistic and special relativistic simulations gets amplified due to the nature of the g 
factor in eqn. (57) which includes the gravitational redshift factor (1 — 2 M/R) 1 ^ 2 and light bending effects. 


4.3. Fully relativistic funnel model 


We construct a general relativistic (Schwarzschild geometry) jet model which consists of a relativistic 
blob in a funnel geometry as shown in Fig. [5] The shape of the magnetic surface can be determined by solv¬ 
ing the relativistic Grad-Shafranov equation (e.g. Fendt & Memola |2001 1 . This is expected to yield stable, 
axi-synmietric magnetic field configurations. Here, we make an approximation of this shape and impose the 
condition of angular momentum conservation along this surface. Here too, we consider the general rela¬ 
tivistic effects which include the time delay due to light bending, the Doppler and gravitational redshift and 
the treatment of the poloidal velocity as a general function of the geometrical and kinematical parameters as 
opposed to a constant that is assumed in the special relativistic case. The kinematical prescription includes 
expressions for the cylindrical radius w and the associated velocity w in terms of 2 and z. The cylindrical 
distance w(z) is given by 

m(z) = wo(l + k(l — e~ z / Zf )), (80) 


where k = (tUf — vjq)/wq, zuj is the cylindrical distance between the normal axis and the source position 
at the location where the funnel transitions into a cylinder. If we use Wf = qwL, 


kf = (Wf - vj 0 )/wq = q/f - 1 . 


( 81 ) 
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Fig. 8.— Simulated light curve for u p = 3 ,xa = 0.9, i = 0.05°, M. = 5 x 10 7 M Q , a = 0.8, wq = 10 tzjl 
and T = 1. The light curve shows the quasi-periodic oscillation expected from a general BL Lac object. A 
general relativistic light curve from orbital features in a cone geometry (blue curve) is compared with the 
previously calculated special relativistic light curve from orbital features in a cone geometry (red curve). The 
increased amplitude ratio and the phase lag change can be attributed to the integrals over time of intensity 
and phase amplifying the small but early differences in phase. Even a small phase difference between the 
general relativistic and special relativistic simulations gets amplified due to the nature of the g factor in eqn. 


(57 1 which includes the gravitational redshift factor (1 — 2 M //() 1/2 and light bending effects. See the online 


article for the color version. 


A constraint onzj can also be obtained based on the vertical distance z at the region where w = Wf. If the 
jet half opening angle is 6 o, tanf?o = (ro/ — wq)/zj. Then 


Zf = (q- f)w L / tan 0 O . 
The velocity associated with w(z) is given by 


. VU()k —zlzt 
■no = z -e 1 f 

z f 

The radial distance R(z) is given by 

R(z) = (w 2 (z) + z 2 ) 1 ' 2 = (tugj 1 + k( 1 - e _2 / Z/ ) 2 + z 2 ) 1 ' 2 

and the velocity associated with R(z) is given by 

R = -(mw + zz) = -(^-(l + k(l- e- z l z f)e~ z l z i + z 
R R \ Zf 


(82) 


(83) 


(84) 


(85) 


As before, the polar angle 9(z) is given by 


sin0(^) = w{z) / R{z) 


( 86 ) 
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and the velocity associated with 9{z) is given by 


M = 1 U - If = + !)-(!+ k)). 


(87) 


The conservation of angular momentum gives an azimuthal velocity 0 = fi(z) = joo /ro 2 . The condi¬ 
tion u ■ u = — 1 in case of the Schwarzschild metric gives 


u* = 


(1 - 2 M/R) - 


j 2 

Joo 

w 2 c 2 


R 2 c 2 


Wq k , 

z f 


(1 + k{ 1 - e~ z l z s))e~ z l z i + z ) + wl{ke~ z / z f{z/z f + !)-(! + k)f 


( 88 ) 

-1/2 


Using the general expression for u from eqn. (74 1 , we obtain an equation for z given by 


z = cRl(l- 2 M/R) - 


i 2 

Joo 


(1 - 2M/w 0 ) 2 1 


1/2 


W 2 C 2 (1-2 M/R) 2 7 % t .J 

{ 2 h, \ 2 

(^(1 + k( 1 - e~ z / z f))e~ z / z f +zj + wl{ke~ z / z f{z/z f + !)-(! + k)f 


- 1/2 


(89) 


In the above equation, w = w(z) from eqn. (80 1 and R. = R(z) from eqn. (84). The above expression 
can be integrated to obtain z = z(t) as 


^ dz (^(! + *(1 - e~ zlz f))e~ z / z f + z^j + w 2 (ke~ z ' z f (z/z f + !)-(! + k )) : 




1/2 


R 


\ 


(1 - 2 M/R) - 


jL (1 - 2 M/wqY 1 


(1 - 2M/R) 2 


f 

Jo 


= c dt. 




(90) 


If we use k from eqn. (81) and Zf from eqn. (82», we can obtain z = z(t) which can then be used to 
evaluate i = z(z(t)), w = w(z(t)), w = w(z(t)). The magnitude of the source velocity is given by 


|i,| = (X, . x.7 2 = + i 2 ) 1/2 = (i 2 + 1) + jl/w 2 \ 


1/2 


(91) 


Once we obtain 2 = z(t), this can be used to evaluate z = z(z(t)), w = w(z(t )), w = w(z(t)). The 
velocity i = z(z(t)) and all other terms in the above expression can be cast in terms of z(t). In the current 


model, we consider the light bending effect in which case the direction vector ko is obtained from eqn. (37 1 
and the angle between the velocity vector x s and the initial direction vector ko given by cos £ is given by 


the general expression in eqn. (55 i. 
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The effective redshift factor g(t) is then given by the general expression in eqn. (57 ). We can use cos£ 
from eqn. (55 I, 7 j e t from eqn. (74 1, (3j et from eqn. (36> in the equation for g to obtain g = git). The light 
curve is given by F(t 0 b s ) = g 2+r {t 0 b s ), expressing t in terms of t n f )S . In this case too, t„i ls tends to increase 
systematically with t em analogous to the full relativistic cone case. Though, the change is slightly gradual. 
This difference in timescale is of the order of ~ 3 days by the end of the simulation. This difference is 
plotted in Fig. [9j clearly the use of the fully relativistic model impacts the phase significantly. 



Emission time r em (days) 


Fig. 9.— Curves showing t 0 b s versus t em for the fully relativistic funnel model (black) and the special 
relativistic cone model (red, dashed). The t 0 b s in the former increases systematically and reaches a difference 
of ~ 3 days by the end of the simulation with respect to the latter. See the online article for the color version. 


To illustrate the model developed in lj3] and the comparison with the special relativistic cone case in 
§4.1[ we simulate a light curve for a general BL Lac object using the expressions presented in this section. 
Using once again u p = 3 as an initial condition, i = 0.05°, M, = 5 x 10 7 Mq, a = 0.8 and mo = 10m l, the 
light curve is plotted in Fig. 10 The parameter u p = 3 is chosen in order that "fjet.i at large m ~ 3.2 which 
is within the range predicted by the effects of radiation pressure and drag in { 2.1 and is in the same range 
as that considered in Camenzind & Krockenberger ( [1992 >. A comparison between the light curves in the 
general relativistic cone case and the special relativistic case indicates an increased amplitude in the former 
in the funnel geometry as well; there is also a systematically increasing phase lag compared to the latter case 
as expected due to the effect of time delay caused by light bending. The ratio of the amplitudes gives a 12 % 
increase in the maximum amplitude due to the GR boost factor in g justifying the use of a fully relativistic 
model. Thus, the beaming effect is observed for the first two cycles for the general relativistic funnel case 
where there is an increased amplitude which then settles down. The increased amplitude ratio and the phase 
lag change here too can be attributed to the integrals over time of intensity and phase amplifying the small 
but early differences in phase. Even a small phase difference between the general relativistic and special 
relativistic simulations gets amplified due to the nature of the g factor in eqn. (57) which includes the 
gravitational redshift factor (1 — 23///s') 1 /2 and light bending effects. 


A comparison is then done between the simulated light curve in funnel and cone geometries for the 












-29- 



Fig. 10.— Simulated light curve for u p = 3 ,xa = 0.9, i = 0.05°, AT. = 5 x 1O 7 M 0 , a = 0.8, wq = 10 zul 
and T = 1. The light curve shows the quasi-periodic oscillation expected from a general BL Lac object. A 
general relativistic light curve from orbital features in a funnel geometry (black curve) is compared with the 
previously calculated special relativistic light curve from orbital features in a cone geometry (red curve). The 
increased amplitude ratio and the phase lag change can be attributed to the integrals over time of intensity 
and phase amplifying the small but early differences in phase. Even a small phase difference between the 
general relativistic and special relativistic simulations gets amplified due to the nature of the g factor in eqn. 


(57 1 which includes the gravitational redshift factor (1 — 2AT//() 1/2 and light bending effects. See the online 


article for the color version. 


general relativistic formalism. For the parameter values jjet,i = 4, k = 2 (q/f = 3), 9 0 = 0.1°, i = 5°, 
AT. = 5 x 10 7 AT.,, a = 0 and zuq = 10 vjl, the light curves are plotted in Fig. 11 The phase of the 


funnel geometry curve systematically lags the phase of the cone geometry curve indicating a dominance by 
higher frequencies. This phase shift is seen to increase with the simulation time. The ratio of the maximum 
amplitude of the funnel geometry curve to that of the maximum amplitude of the cone geometry curve is 
9 %. Because of the slower expansion of the funnel compared to the cone there are more orbits at higher 
frequencies (due to conservation of angular momentum) at a location where the GR boost factor is more 
effective. 


4.4. Funnel model simulations and discussion 


We perform two sets of simulations to address the expected light curve and timing information for an 
observer viewing the AGN at varying inclination angles i, black hole mass AT. and the initial launch bulk 
Lorentz factor "fjet.i- The choices for i are based on typical values expected for the inclination angle from 


observational studies (e.g. Pushkarev et al. 2009 

). A time series analysis of the light curve is conducted 

using the Fourier periodogram analysis (e.g. 

Mohan et al. 2014, 

Mohan & Mangalam 2014) and the wavelet 
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Observer time / ohs (days) 


Fig. 11.— Simulated light curve for 7 j e t,i = 4, k = 2, xa = 0.9, i = 0.05°, M. = 5 x 1O 7 M 0 , a = 0 
and tz7(j = lOru/,. A comparison is made between the funnel geometry light curve (black curve) and the 
cone geometry light curve (blue curve). Both light curves show a quasi-periodic oscillation expected from 
a general BL Lac object. The phase of the funnel geometry curve systematically lags the phase of the cone 
geometry curve. Also, the funnel geometry curve is dominated by higher frequencies in the initial portion. 
See the online article for the color version. 


analysis (e.g. Torrence & Compo 1998; Moh an et al.|2011[ ). The Fourier periodogram of the light curve is the 
normalized Fourier power spectrum, evaluated at the frequencies / = j/(N At) where j = 1,2, ...(N/2— 1) 
(upto and excluding the Nyquist frequency) and is ht with a power law shape, P(f) = Af m , assumed to 
be the shape that best describes the underlying power spectral density (PSD). The wavelet analysis gives 
the QPO and its evolution including the phases during which it is present in the light curve and hence, 
the number of cycles it is present for. In the wavelet analysis, any periodicity inside a cone of influence 
(triangular region) can be trusted. Features outside this region could be subject to systematic effects due to 
the wavelet method. From these simulations, we aim to address questions such as the typical QPO expected 
from orbital processes in the jet, its dependence on M., the beaming effect, the conditions for its sustenance 
and the typical range of PSD slopes expected. These are motivated from observational studies at optical and 


radio wavelengths which indicate QPOs with timescales of less than a day to a few tens of days (e.g. Gupta 
et al.| 2012 ] ). 


We first simulate the expected light curve from an emitting source in helical motion along the funnel 
shaped magnetic surface launched from a single ring with cylindrical radius w = lOzu^. The other quanti¬ 
ties which are fixed include k = 2, xa = 0.9, 9q = 0.1° and a = 1. The simulations are carried out for i = 
3°, 6 °, 9°, 12° and 15° with M, = (0.5, 5) x 1O S M 0 and 7 j etj i = 2 (mildly relativistic), 4 (relativistic) and 
10 (highly relativistic). The choice of 7 j e t,i = 2,4 is consistent with the range predicted by the effects of 
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radiation pressure and drag in §2.1| The choice of 7 , = 10 was made in order to study simulated light curves 
for jets from black holes with higher masses which includes blazars. The results of the timing analysis are 
summarized in Tableland some interesting cases are plotted in Fig. 

With increase in the inclination angle i, the maximum amplitude in the simulated light curve tends 
to decrease. This is due to the beaming effect as the projection of the emitting source velocity onto the 
emission direction (/3j e t ■ k) continues to grow smaller in the direction towards the observer. This leads to 
a less pronounced maximum amplitude which continues to decrease with increase in i. For a given M,, 
the maximum amplitude tends to increase with "( 3e t,i- This is due to the stronger beaming effect towards 
the observer line of sight by the emitting source for larger 7 j e t,i- The PSD slope ranges between —1.54 
and —3.43. The slopes steeper than —2.5 are likely to be the result of a poor fit. Hence, we calculate a 
median slope of —2.00 ± 0.28. QPOs peaked between 20.7 - 130.7 days are inferred from the wavelet 
and PSD analysis. For a given M», the QPO peak shifts to higher values with increasing i. This is due to 
the orientation away from the observer’s line of sight for the chosen set of parameters leading to the larger 
timescale QPO at larger r. 

We then simulate the expected light curve from multiple emitting sources in helical motion along the 
funnel shaped magnetic surface launched from rings with cylindrical radii w = 5, 6 , 7 and 8 wl- Fixed 
parameters in this simulation include k = 2 , xa = 0.9, do = 0.1° and a = 1. The simulations are earned 
out for i = 3°, 6 °, 9°, 12° and 15° with M, = (0.5, 5) x 10 8 M Q and 7 j ett i = 2 — 10 as in the previous case. 
The results of the timing analysis are summarized in Table [2] and some interesting cases are plotted in Fig. 

M 

In these simulations too, we observe similar trends as were present in the previous case. With increase 
in i, the maximum amplitude in the simulated light curve tends to decrease. For a given M„ the maximum 
amplitude tends to increase with 7 j e t,i- The PSD slope in these simulations is more well constrained and 
ranges between —1.09 and —3.39. The slopes steeper than —2.5 are likely to be the result of a poor fit. 
Hence, we calculate a median slope of —2.02 ± 0.34. Significant QPOs ranging between 1.37 - 33.49 days 
are detected in these simulations. The QPO for high 7 j e t,i tends to be lower than that for lower 7 j e t,i- This 
is likely as for high 7 j e t, the beamed portion occurs in the inner jet close to the black hole where the orbital 
frequency is higher due to gravitational redshift. Thus, the associated timescales are then expected to be 
lower. Thus, the power spectrum is dominated by higher frequencies which implies that the PSD would tend 
to flatten. This trend is also seen in the above multiple emitting ring simulations for M, = 5 x 10 7 M Q . For 
the case of M, = 5 x 10 8 Mq, the opposite trend is seen. For a higher M„ there are more cycles during 
beaming due to favourable orientation (at higher frequencies); hence, the PSD slope tends to steepen. 


12 


5. Summary & Conclusions 

We summarize our results. The special relativistic cone model (§4.1[) is the same as that proposed by 


Camenzind & Krockenberger|( |l992| ). The fully relativistic models developed in this work include the cone 


model (jj 8 ]l and the funnel model (jjTOj). The main results include discussions on the novel aspects of the 
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Fig. 12.— Results for light curves from a single emitting ring of = 10 vjl- The plots in the left column 
are the simulated light curves, the plots in the middle column are their wavelet analysis and the plots in the 
right column are their PSD. The PSD is fit with a power law with a slope m. The line above the best fit 
indicates a 99% significance on any inferred QPO period. Top row: simulation for i = 3°, M, = 5 x 1O 7 M 0 , 
7 j et = 2; a QPO of 124 days is inferred from the timing analyses; a PSD slope of —1.96 is inferred. Middle 
row: simulation for i = 9°, M, = 5 x 10'M 0 , 7 ^ = 2; a QPO of 126.2 days is inferred from the timing 
analyses; a PSD slope of —2.78 is inferred. Bottom row: simulation for i = 15°, M, = 5 x 1O 7 M 0 , 
7 j et = 4; a QPO of 97 days is inferred from the timing analyses; a PSD slope of —2.50 is inferred. See the 
on line article for the color version. 

fully relativistic funnel model with respect to the special relativistic cone model. 


1. We have considered a possible mechanism of jet based variability from various sources including 
BL Lac objects and blazars (e.g. Camenzind & Krocken berger 1992; Steffen et al. 19951, quasars 
(e.g. Kudryavtseva et al. 20lTj), binary black holes (e.g. |Iguchi et al.|2010j ) and X-ray binaries (e.g. 
Fukumura & Kazanas|2008| Fukumura et al.|2009). We extended the light-house model proposed in 


Camenzind & Krockenberger (1992) which was applied to simulate optical light curves from BL Lac 


objects and quasars. 
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Fig. 13.— Results for light curves from multiple emitting rings of voq = 5 wl to 8m l- The plots in the left 
column are the simulated combined light curves. The plots in the middle column are their wavelet analysis 
and the plots in the right column are their PSD. The PSD is fit with a power law with a slope m. Top row: 
simulation for i = 3°, M, = 5 x 1 O 8 M 0 , 7 j et = 10; a QPO of 1.37 days and a PSD slope of —1.77 are 
inferred from the timing analyses. Middle row: simulation for i = 6 °, M, = 5 x 10 8 Mq, 7 jet = 4; QPOs 
of 4.37 days and 8.74 days (harmonics) and a PSD slope of —2.12 are inferred from the timing analyses. 
Bottom row: simulation for i = 15°, M, = 5 x 1O 7 M 0 , 7 j et = 2; a QPO of 5.45 days and a PSD slope of 
—2.94 are inferred from the timing analyses. See the online article for the color version. 

2. We studied the effects of radiation pressure and drag and derived the saturation Lorentz factors that are 
achievable. The Lorentz factors 7 range between 2 — 7 in the purely radial outflow simulations. For 
a poloidal flow {(5^ = 0), we obtain 7 in the range 1.1 — 26.3, greater than that obtained previously. 
Though, the drag force acting on the non-radial component of the blob velocity plays an important 
role in rapidly decreasing 7 ~ 2 thus stabilizing it at very small x, indicating that the outflow is stable 
even in the innermost regions. Thus, the 7 j e t,i chosen for the simulations in the relativistic cone and 
funnel models in the subsequent section can be obtained at larger distances with small to moderately 
relativistic emitting blob velocities. 
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3. We constructed a kinematic model of orbital blobs along helical trajectories on a magnetic surface 
approximating the expanding jet with foot points on the accretion disk in the vicinity of a black hole 


using a special relativistic calculation of the g factor for cone geometry similar to (Camenzind & 
Krockenberge r|1992 1 in f 4.1 and a fully relativistic formulation in cone geometry ( ; J 8 ]) and the same 
in a more realistic funnel geometry (^T0]>. The g factor was calculated in Schwarzschild geometry and 
the following GR effects were included: gravitational and Doppler shifts, aberration and a prescription 
for time delay due to orientation and light bending. 


4. By using the periodogram and wavelet based time series analysis techniques ( ; |4.4[ ), we seek to dis¬ 
tinguish amongst the various flavours of the generic jet variability model which here include the cone 
and funnel models. The time series analysis of simulated light curves yields properties including ape¬ 
riodic variability, its timescales and the emission region; the power spectral density slope and its range 
which can be compared with observations in optical/UV and X-ray; it is also found that the QPO, its 
evolution and stability can be described by the orbital motion of blobs in the jet. 


5. The use of fully relativistic models (given in (| 8 ]and ' |T()| ) was justified by the resulting amplitude 
increase (by about 12 %) due to the GR boosts and a systematically increasing phase lag in the 
simulated light curve for a general BL Lac object from the fully relativistic funnel when compared to 
the special relativistic cone which occurs at small k and larger inclination angles i. The phase lag has 
been explained by the angular momentum conservation in combination with the gradual increase of 
the orbital radius. Our fully relativistic formulation reduces to the special relativistic formulation in 
the limit of large R and when ah the above GR effects are not considered; this is presented in J4] 


6 . There is thus a necessity to account for all GR effects as used in the present model for the correct 
description of the physical effects on emitted radiation from these jetted sources and to use more 
realistic jet geometries. It can be thus applied to both timing studies of jet variability as well as to map 
trajectories of radio blobs or blobs in the inner jet and compared with observations to help identify the 
region of emission. The application of the model to radio data is being planned. 


7. Two sets of simulations, one for blobs emanating from a single and another from multiple rings were 
carried out to span a range of light curves. Their - timing properties were studied for the maximum 
amplitudes, PSD slope and QPO if present. Two main trends are observed in the single emitting 
ring simulations. With increase in the inclination angle i, the maximum amplitude in the simulated 
light curve tends to decrease, as the beaming effect in the direction of the observer continues to grow 
smaller. The other trend is that for a given M„ the maximum amplitude tends to increase with 7 
due to the stronger beaming effect towards the observer line of sight by the emitting source for larger 
"fjet.i- We calculate a median slope of —2.00 ± 0.28 similar to the observed PSD slopes of blazars 
in the Optical and X-ray translated to the source frame. QPOs peaked between 20.7 - 130.7 days are 
inferred from the wavelet and PSD analysis. The second set of simulations were carried out for an blob 
launched from multiple emitting rings wq = 5 — 8m l- In these multiple emitting rings simulations 
too, the same above two trends are observed. We calculate a median slope of —2.02±0.34. Significant 
QPOs ranging between 1.37 - 33.49 days are detected in these simulations. In the multiple emitting 
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rings simulations, as we were also able to identify QPOs as a wider choice of orientations result from 
a larger set of initial conditions; same trends were observed in the relationship between the QPO 
timescale, M, and 7 j et . The QPO timescale tends to reduce with increase in 7 j et due to GR boost 
which tends to flatten the PSD slope. Aperiodic variability when considering multiple emitting rings 
is a natural consequence that is seen in our results. 

A natural power law shaped PSD with a typical slope of ~ — 2 along with a weak to strong QPO 
ranging between 1.37 - 130.7 days emerges from the simulations considering single and multiple emitting 
rings. This shape comes naturally considering only geometrical parameters, even before any consideration 
of physical models of the orbital instabilities. The detailed magnetic structure can be constructed by solving 
the relativistic Grad-Shafranov equation. The dynamical timescale due to jet based orbital processes is 
expected to be from a few minutes to tens of days depending on the black hole mass. The general relativistic 
funnel model for variability can thus be applied to radio, optical and X-ray emission from various types of 
jetted sources including radio loud AGN such as BL Lac objects, quasars, binary black holes, X-ray binaries 
where active accretion results in the jet phenomenon such as in micro-quasars and other compact sources 
such as accreting and hence active neutron stars. 
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A. Relationship between the effective redshift factor g and the Doppler factor D 

The tetrads used by a local static observer in the Schwarzschild metric are given by, 


ef t) = (1 - 2Af/ J R)~ 1 / 2 (l, 0,0,0) (Al) 

ef r) = (l-2M/ J R) 1 /2( 0 ,1,0,0) 

^ = (1/^(0,0,1,0) 

e «>) = (!/-R sin 0 ) (0,0,0,1). 


The directional vector rp a) in the general relativistic case with light bending is given by (Pechacek et al. 


2005), 


n ^ = {n^ t \n^ r \n^\n^ >S> ) = ( 1, cos a,-- cosi, — --- sin^sinl ) . (A2) 

\ sm ip sm ip J 

If p is the four-momentum of the emitted light ray, its covariant components are given 

Pn=Ptil,Pr/pt,Pe/Pt,P<i>/pt), (A3) 

the objective being to calculate the components pp in terms of the emission angles and the geometric quan¬ 
tities. This can be done using the relationship between the directional vector n ta> and given by 


B W = 


P» e \a) 


(A4) 


Using the above equation and the tetrads from eqn. m we obtain 

Pr/pt = cos a(l — 2M/R)~ 1 (A5) 

Pe/pt = — S ^ na cosiR(l — 2M/R)- 1 / 2 
am ip 

Pch/pt = — Sma sin 6 sin iR sin 9(1 — 2M /i ?) -1 / 2 
am ip 


The four velocity of the emitting source is given by 




(A 6 ) 


where u 1 = 7 j et is obtained from eqn. (35 1 . The effective redshift factor g is the ratio of observed to emitted 
energy of the photon and is given by, 


-^observed _ P(®®) ' u(oo) _ 1 

-Remitted P ’ U ' 


(A7) 
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Using eqns. (A5 I and (A 61 in the above equation, the effective redshift factor is given by, 


(1 - 2M/R) 1 / 2 
u l 


sin a 
sin if; 


R/3g cos i — 


((\ — 2M/R) 1 / 2 — cos a- - 

V (1 — 2M/R) 1 / 2 

sin Oi . , . . 

-- R sm tip a, sin sin i 

sin ip 



(AS) 


The g expression can be written in terms of (3j e t cos£ with d je t obtained from eqn. (361 already 
including a correction for the frame changing from an instantaneous rest frame to a local static observer 


frame and cos £ obtained from eqn. (55 l to give the expression for g in eqn. (57 1 . Thus 


9 = 


(1 - 2 M/R) 1 / 2 
ljet {1 - Pjet COS O’ 


(A9) 


As p M can also be expressed as 

Vll = e(l, b(l/b 2 - 1/R 2 (1 - 2M/R))- 1 / 2 ,p e /e,b) 1 (A10) 


where e is the conserved photon energy, the direction between the emitted ray and the radial direction for 
equatorial orbits (9 = tt/2) is given by 


P" e W _(Vfl)_ 

Pntf r) (1 - 2M/R))- 1 / 2 b{l/b 2 - l/R 2 (l - 2M/R))~ 1 / 2 ' 


Solving the above equation for sin a, we obtain eqn. (42 1 . 


(All) 


REFERENCES 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011, ApJ, 726, 43 
Abramowicz, M. A., Ellis, G. F. R., & Lanza, A. 1990, ApJ, 361, 470 
Beloborodov, A. M. 2002, ApJ, 566, L85 

Brandenburg, A., Nordlund, A., Stein, R. F., & Torkelsson, U. 1995, ApJ, 446, 741 

Camenzind, M., & Krockenberger, M. 1992, A&A, 255, 59 

Doeleman, S. S., Fish, V. L., Schenck, D. E., et al. 2012, Science, 338, 355 

Fendt, C., & Memola, E. 2001, A&A, 365, 631 

Fukumura, K., & Kazanas, D. 2008, ApJ, 679, 1413 















-38- 


Fukumura, K., Kazanas, D., & Stephenson, G. 2009, ApJ, 695, 1199 

Gaur, H., Gupta, A. C., Lachowicz, P., & Wiita, P. J. 2010, ApJ, 718, 279 

Gupta, A. C., Fan, J. H„ Bai, J. M„ & Wagner, S. J. 2008, AJ, 135, 1384 

Gupta, A. C., Krichbaum, T. P„ Wiita, P. J., et al. 2012, MNRAS, 425, 1357 

Hagen-Thorn, V. A., Larionov, V. M., Jorstad, S. G., et al. 2008, ApJ, 672, 40 

Horak, J., & Karas, V. 2006, PASJ, 58, 203 

Iguchi, S., Okuda, T., & Sudou, H. 2010, ApJ, 724, L166 

Kudryavtseva, N. A., Britzen, S., Witzel, A., et al. 2011, A&A, 526, A51 

Lachowicz, R, Gupta, A. C., Gaur, FL, & Wiita, P. J. 2009, A&A, 506, L17 

Liu, H„ Wang, J., Mao, Y., & Wei, J. 2010, ApJ, 715, LI 13 

Mangalam, A., & Mohan, P. 2014, Journal of Astrophysics and Astronomy, 35, 431 

Mangalam, A. V., & Wiita, P. J. 1993, ApJ, 406, 420 

Markoff, S., Nowak, M. A., & Wilms, J. 2005, ApJ, 635, 1203 

Marscher, A. P, & Gear, W. K. 1985, ApJ, 298, 114 

Marscher, A. P, Jorstad, S. G., D’Arcangelo, F. D., et al. 2008, Nature, 452, 966 
Mohan, P. 2014, PhD Thesis, Indian Institute of Astrophysics, Bangalore. 

Mohan, P, & Mangalam, A. 2014, ApJ, 791, 74 

Mohan, P, Mangalam, A., Chand, H., & Gupta, A. C. 2011, Journal of Astrophysics and Astronomy, 32, 
117 

Mohan, P, Mangalam, A., & Chattopadhyay, S. 2014, Journal of Astrophysics and Astronomy, 35, 397 

Novikov, I. D., & Thorne, K. S. 1973, in Black Holes (Les Astres Occlus), ed. C. Dewitt & B. S. Dewitt, 
343-450 

Pechacek, T., Dovciak, M., Karas, V., & Matt, G. 2005, A&A, 441, 855 
Poutanen, J., & Beloborodov, A. M. 2006, MNRAS, 373, 836 

Pushkarev, A. B., Kovalev, Y. Y., Lister, M. L., & Savolainen, T. 2009, A&A, 507, L33 
Qian, S. J., Quirrenbach, A., Witzel, A., et al. 1991, A&A, 241, 15 



-39- 


Rieger, F. M. 2004, ApJ, 615, L5 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Spruit, H. C. 1996, ArXiv Astrophysics e-prints, arXiv:astro-ph/9602022 

Steffen, W„ Zensus, J. A., Krichbaum, T. P„ Witzel, A., & Qian, S. J. 1995, A&A, 302, 335 

Torrence, C., & Compo, G. R 1998, Bulletin of the American Meteorological Society, 79, 61 

Vokrouhlicky, D., & Karas, V. 1991, A&A, 252, 835 

Wagner, S. J., Witzel, A., Heidt, J., et al. 1996, AJ, 111, 2187 

Wiita, R J. 2006, in Astronomical Society of the Pacific Conference Series, Vol. 350, Blazar Variability 
Workshop II: Entering the GLAST Era, ed. H. R. Miller, K. Marshall, J. R. Webb, & M. F. Aller, 
183 

Zhang, X.-H., & Bao, G. 1991, A&A, 246, 21 


This preprint was prepared with the AAS IATpX macros v5.2. 



-40- 


Inclination 

Angle 

*(°) 

M. 

'ijetj, 

Max. 

Amplitude 

A 

PSD 

Slope 

m 

QPO 

(days) 

10 8 M 0 

3 

0.5 

2 

54.48 

-1.96 

124.0 W 



4 

506.49 

-1.97 

89.8 W 



10 

8090.84 

-2.26 

79.7 W 


5 

2 

54.24 

-1.54 

68.2 



4 

492.43 

-1.65 

18.6 



10 

6949.31 

-1.82 


6 

0.5 

2 

54.05 

-2.42 

124.8 W 



4 

483.05 

-2.18 

90.8 W 



10 

5780.68 

-3.00 

80.8 W 


5 

2 

54.06 

-1.64 

70.1 



4 

483.10 

-1.73 

20.7 



10 

5780.68 

-2.01 


9 

0.5 

2 

51.50 

-1.70 

126.2 W 



4 

372.24 

-2.30 

92.4 W 



10 

1540.67 

-3.16 

82.5 W 


5 

2 

51.51 

-1.70 

73.2 



4 

372.25 

-1.77 

24.3 



10 

1544.25 

-2.03 


12 

0.5 

2 

46.60 

-1.73 

128.2 W 



4 

238.24 

-2.42 

94.8 



10 

328.10 

-3.13 

85.0 W 


5 

2 

46.60 

-1.73 

77.5 



4 

238.24 

-1.80 

29.0 



10 

328.56 

-1.98 

- 

15 

0.5 

2 

40.24 

-3.43 

130.7 W 



4 

136.38 

-2.50 

97 W 



10 

81.36 

-3.07 

88.1 W 


5 

2 

40.25 

-1.75 

82.9 



4 

136.38 

-1.82 

35.1 



10 

81.44 

-1.94 



Table 1: Results for light curves from blobs launched from a single ring of cylindrical radius zuq = 10z&l. 
The parameters used which arc not reported in the above table include x a = 0.9, 6$ = 0.1° and a = 1. 
The slopes steeper than —2.5 arc likely to be the result of a poor fit. Hence, we calculate a median slope of 
—2.00 ± 0.28 and the QPO estimates range between 20.7 - 130.7 days. A W next to a QPO peak indicates 
its measurement using the wavelet analysis only. If there is nothing present next to the inferred QPO, the 
measurement is the mean of the wavelet and PSD estimates. 
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Inclination 

Angle 

n°) 

M. 

r Ijet,i 

Max. 

Amplitude 

A 

PSD 

Slope 

m 

QPO 

(days) 

10 8 M 0 

3 

0.5 

2 

54.35 

-2.01 

- 



4 

494.97 

-1.67 

- 



10 

5955.14 

-1.70 

- 


5 

2 

47.94 

-1.72 

27.62 



4 

327.01 

-1.90 

7.89 



10 

1909.06 

-1.77 

1.37 

6 

0.5 

2 

52.28 

-2.47 

- 



4 

424.37 

-1.73 

- 



10 

4155.13 

-1.96 

0.98 


5 

2 

47.43 

-1.95 

28.34 



4 

338.99 

-2.02 

4.37, 8.74 



10 

3553.40 

-2.36 

1.81 

9 

0.5 

2 

48.46 

-2.81 

- 



4 

322.86 

-1.81 

- 



10 

2778.43 

-1.31 

- 


5 

2 

45.01 

-2.12 

14.80, 29.59 



4 

298.02 

-2.18 

5.06, 10.13 



10 

2778.43 

-3.18 

- 

12 

0.5 

2 

44.00 

-2.99 

- 



4 

217.66 

-2.00 

- 



10 

962.85 

-1.09 

- 


5 

2 

40.96 

-2.39 

15.65, 31.31 



4 

218.57 

-2.34 

6.03 



10 

962.85 

-3.39 

- 

15 

0.5 

2 

38.23 

-2.94 

5.45 



4 

135.59 

-1.97 

- 



10 

207.60 

-1.23 

- 


5 

2 

35.78 

-2.61 

16.74, 33.49 



4 

137.58 

-2.46 

2.89, 7.22 



10 

207.60 

-3.24 

1.89 


Table 2: Results for light curves from blobs launched from multiple rings of cylindrical radii vjq = 5 — 8m l. 
The parameters used which arc not reported in the above table include x a = 0.9, 6q = 0.1° and a = 1. 
We calculate a median slope of —2.02 ± 0.34 and the QPO estimates range between 1.37 - 33.49 days. All 
QPOs inferred above are the mean of the wavelet and PSD estimates. 




